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Abstract 


Various implementations of the direct simulation Monte Carlo 
(DSMC) method exist in academia, government and industry >. By 
comparing implementations, deficiencies and merits of each can be 
discovered. This document reports comparisons between DSMC Analysis 
Code (DAC) and “MONACO”. DAC is NASA’s standard DSMC 
production code and MONACO is a research DSMC code developed in 
academia. These codes have various differences; in particular , they 
employ distinct computational grid definitions. On one hand, DAC 
employs a wall surface grid comprising unstructured triangular cells in 
conjunction with a volume grid comprising three-space Cartesian cells. 
On the other hand, MONACO employs a volume grid generally 
comprising unstructured cells and defines a wall surface by one of its 
surface boundaries. In this study, DAC and MONACO are compared by 
having each simulate a blunted flat plate wind tunnel test, using an 
identical volume mesh. Simulation expense and DSMC metrics are 
compared. In addition, flow results are compared with available 
laboratory data. Overall, this study revealed that both codes, excluding 
grid adaptation, performed similarly. For parallel processing, DAC was 
generally more efficient. As expected, code accuracy was mainly 
dependent on physical models employed. 


Nomenclature 


D total drag, N 

Dpress drag due to pressure, N 

D S hear drag due to shear, N 

Err error 

L flat plate length 

m-t serial processing memory 

rriproc parallel memory of processor proc 

n ra free stream number density, m -3 

N p number of simulation molecules 

N proc number of parallel processors 

S pm parallel processing memory spread 

t-L serial processing real time, s 

t parallel parallel processing real time, s 

T'ref collision model reference 

temperature K 

Tj ot reference temperature for rotational 

energy exchange, K 


x, y,z 
Z r 


& diameter 

A, As 

£ rot 

trot 

dps 


X 

co 


isothermal wall temperature, K 
free stream temperature, K 
spatial coordinates, m 
rotational collision number 
maximum rotational collision 
number 

collision diameter, m 

change, characteristic cell size 

molecule rotational energy, / 

rotational degrees of freedom 

parallel processing speed up 

efficiency 

mean free path, m 

free stream mean free path, m 

collision time, s 

reference collision time, s 

viscosity index 



1 Introduction 


Various implementations of the direct simulation Monte Carlo (DSMC) method exist in academia, 
government and industry. By comparing implementations, deficiencies and merits of each can be 
discovered. This document reports comparisons between DSMC Analysis Code (DAC) and “MONACO”. 
DAC is NASA’s standard DSMC production code and MONACO is a research DSMC code developed in 
academia. These codes have various differences; in particular, they employ distinct computational grid 
definitions. On one hand, DAC employs a wall surface grid comprising unstructured triangular cells in 
conjunction with a volume grid comprising three-space Cartesian cells. On the other hand, MONACO 
employs a volume grid generally comprising unstructured cells and defines a wall surface by one of its 
surface boundaries. In this study, DAC and MONACO are compared by having each simulate a blunted 
flat plate wind tunnel test, using an identical volume mesh. 

This document begins by providing a brief overview of DAC and MONACO. Then, a description on 
their setup for comparisons is given. Next, a description of a flat plate wind tunnel test is provided. Then, 
boundary conditions and collision model configurations are outlined. Subsequently, computational mesh 
generation and configuration is presented, including a grid convergence study. 

Following these descriptions, comparisons of resulting flat plate simulations are presented. First, serial 
and parallel simulation time and memory are compared, along with code function profiles of serial 
simulations. Second, the quality of the simulations is examined through DSMC metrics. Third, flow field 
contours of macroscopic flow properties are discussed. Finally, surface properties are studied, including 
comparison against existing surface pressure and heat flux data found in literature. 


2 Description of DSMC Codes 

DAC was invented in the mid 1990’s by LeBeau and others. 1 ’ 2 The code is written in FORTRAN 77, 
and uses Cray style pointers to permit dynamic memory allocation. DAC has been employed by NASA 
and other organizations for mission support, aerospace vehicle design analysis, computational 
nonequilibrium gas modeling and in support of laboratory projects. 3 Simulations have involved 
nonequilibrium flow about spacecraft in orbit, 4 ” 6 rarefied hypersonic flow 7 ’ 8 and support of hybrid 
simulation research. 4 ’ 8 For DAC, each computational grid is a composition of a vehicle surface grid and a 
flow-field volume grid. By using two grids, the surface grid does not impose a geometric constraint on the 
surrounding volume grid. The vehicle surface grid, generated independently, is a closed surface 
comprising an unstructured distribution of triangular cells. The flow-field grid is generated automatically 
by DAC’s preprocessing utility called PREDAC. It is nominally bounded by six planes, each coordinate 
axis perpendicular to a pair, and by the vehicle surface grid. Each bounding plane represents an inflow, 
outflow or symmetry boundary. Alternatively, the computational domain can be defined more generally 
via the use of the same type of unstructured triangular grids that are used to define a surface, as the 
boundary conditions for these grids can be defined as surface, inflow or outflow. The initial flow-field 
grid comprises three-space Cartesian cells, labeled “level 1” cells, often distributed uniformly along each 
coordinate direction for initial simulation. Solution results obtained from level 1 cells enable PREDAC to 
refine the grid in regions of higher flow density by dividing level 1 cells into smaller cells, labeled “level 
2” cells. The version of DAC employed for this study is based on 97 distribution 11. 

MONACO is a research code developed in academia. Its name is a geographic reminder of Monte 
Carlo simulation. The code was devised in the mid 1990’s by Dietrich and Boyd. 9 ' 10 The code is written in 
the C programming language. MONACO has been employed by research universities and laboratories for 
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computational nonequilibrium gas modeling and in support of laboratory projects. Simulations have 
involved support of hybrid simulation research, “ nonequilibrium hypersonic flow “ and 
nonequilibrium flow about spacecraft. 23 25 MONACO employs a single body fitted grid with general 
surface topology for its boundaries. External software is used to generate the grid and define inflow, 
outflow, symmetry and wall surface boundaries. Two or three-dimensional grids are applicable with 
structured or unstructured cell distributions. This study employs three-dimensional grids with Cartesian 
cells. Mesh adaptation to simulation results is performed with external software. The version of 
MONACO employed for this study is based on 3.09. 

For this study, both codes employ DSMC algorithms based on Bird’s 1994 algorithms, 26 with random 
pairing of molecules within a cell for collision. Gas molecule collision dynamics follow the variable-hard- 
sphere (VHS) 27 model. Both codes consider rotational energy exchange according to Borgnakke-Farson 
statistical model. 28 However, DAC assumes constant rotational collision number, while MONACO 
employs the variable rotational energy exchange probability model of Boyd. 29 For gas-surface 
interactions, both codes are run with diffuse reflection and full thermal accommodation, against an 
isothermal wall. Gridgcn ’° is employed to generate surface grids for DAC and volume grids for 
MONACO. 


3 Description of Flat Plate Simulations 

Wind Tunnel Model 

Notwithstanding differences in grid processing procedures between DAC and MONACO, flow field 
computational grids exist that are identical in both codes. Excepting solid surface discretization, these 
grids enable grid independent comparisons between the codes. For this study, such comparisons are 
achieved by simulation of rarefied flow about rectangular flat-plate model within a free-jet expansion 
wind-tunnel. The particular geometry and flow conditions are taken from an experiment presented by 
Allegre and others in 1992. 11 This experiment is characterized by a rarefied global Knudsen number of 
0.016 based on the free stream mean-free-path 7 l x = 1.60 mm (computed from free-stream conditions) 
and flat-plate length L = 100 mm. Flat plate thickness is 5 mm. Further details of the wind-tunnel 
experiment are found in references 31-33. 

Flow Conditions 

Identical flow conditions are employed for each code to simulate the experiment. Both codes 
commenced with a vacuum domain. Inflow of pure nitrogen is assumed. The following input parameters 
define the nitrogen gas species: molecular weight MW = 28.01, VHS molecular diameter 8 diameter = 
4.17x10 1 °?ti, 26 viscosity index o> = 0.74 at reference temperature T re / = 273 K, 26 and number of 
rotational degrees of freedom d/ rQt = 2. However, because of distinct rotational relaxation models, DAC 
utilizes constant rotational collision number Z r = 5, while MONACO employs variable rotational 
collision number model with reference temperature for rotational energy exchange T/ ot = 91.5 K and 
limiting rotational collision number Z roo = 18.1. Unless otherwise indicated, simulations are performed 
with molecule rotational energy exchange procedures activated. Nitrogen gas enters the inflow boundary 
uniformly with = 1503 m/s, T oa = 13.32 K and n m = 3.716xl0 20 ?n -3 . Free stream velocity is 
parallel to x axis and longitudinal surfaces of flat plate. The flat-plate boundary condition is an isothermal 
wall temperature T w = 290 K. Gas molecules are assumed to reflect diffusely from the wall according to 
the Maxwellian velocity distribution at T w . 
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Computational Mesh 

The computational flow field grid selected for baseline simulation is a three-dimensional structured 
grid with cubic cells having 0.5 mm sides. The domain is a rectangular region of flow about the flat-plate 
and has two-dimensional symmetry, and is illustrated by Figure 1. The domain dimensions are 180 by 
205 by 1 mm in the x, y and z directions, respectively. Inflow boundaries are defined by left, bottom and 
top faces in Fig. 1. Outflow boundary is defined by the right face. Symmetry boundaries are defined by 
the front and back faces. Flat plate surface boundary is defined by the open interior region with 
dimensions 100 by 5 by 1 mm. in the x, y and z directions, respectively. The origin of the coordinate 
system is at the center of the flat plate left face. With DAC, flow field mesh generated with level 1 cells 
so that MONACO could easily use the same flow field mesh. To keep the focus of the comparisons on the 
core DSMC algorithms, variable particle and time weighting is not employed. 

Flat plate surface discretization is illustrated by Figs. 2 (a) and (b). Here a clear difference between the 
codes is highlighted. DAC’s surface mesh is independent from the flow field mesh and must employ 
triangular cells, whereas, MONAC’s surface mesh is a boundary of the flow field mesh and utilizes 
associated cell faces. 

To determine baseline cell size, simulations were performed with three cubic cell edge sizes: 1, 0.5 
and 0.25 mm. The largest of these is smaller than the free stream mean free path, 1.6 mm, as computed 
from the free stream conditions. The smallest of these is smaller than the smallest mean free path, X ~ 
0.27 mm, throughout the entire domain. While the domain dimensions were kept constant in the x and y 
directions, it varied in the z direction because the number of cells in the z direction was kept constant at 2. 
Grid refinement was monitored by comparing surface property profiles along the flat plate upper surface. 
Surface pressure, x directed shear stress, incident heat flux and reflected heat flux were examined. In 
addition, pressure drag, shear drag and total drag were monitored. 



Figure 1. Flow field mesh: (a) xy view of simulation boundaries, z axis points out of page, surface of 

flat-plate represented by interior rectangle; (b) isometric view enlarged about lower left 
corner 
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(a) DAC surface mesh 



(b) MONACO surface mesh 


Figure 2. Isometric view of flat plate surface meshes near leading edge 


Figure 3 illustrates effects of grid refinement on shear stress and incident heat flux. For brevity, only 
DAC results are shown. MONACO produces similar results. Qualitatively, it is seen that surface 
properties change subtly with selected changes in grid refinement. Percentage changes in results along the 
surface are also plotted. By inspection, the variation of percentage change along the surface is random; 
thus, the most probable percent change is estimated by the mean. For shear stress and reflected heat flux, 
the mean percentage changes along the upper surface are 1.28 and 1.12 %, respectively, in refining from 
1.0 to 0.5 mm cell edges. In refining from 0.5 to 0.25 mm cell edges, these changes are 1.18 and 0.979 
%, respectively. 
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Figure 3. DAC assessment of grid refinement by upper surface properties 
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Table 1 summarizes grid refinement results by DAC and MONACO. Listed are mean percentage 
absolute value of changes of upper surface pressure, x directed shear stress, and incident and reflected 
heat flux. Additionally, percentage changes of drag, including pressure and shear components, are listed. 
Decreasing edge lengths from 1 to 0.5 mm, results in, local surface property, mean percentage changes 
between 1 and 2 %. Decreasing edges from 0.5 to 0.25 mm, results in these changes being less than 1 %, 
except for shear stress. The reason for sensitivity of shear stress is suspected due to selection of only two 
cells along z direction. Now examining total surface drag, it is seen that changes are less than one percent, 
even in refining from 1.0 mm to 0.5 mm cell edges. Because the percentage changes of total surface drag, 
and incident and reflected heat flux, decreases from the first to the second refinements, and because they 
are less than 1 % in the second refinement, the computational mesh with 0.5 mm cubic cell edge lengths 
is selected as baseline for further DSMC code comparisons involving selected simulation case. 


Table 1. Effects of grid refinement on surface properties with the following 

cubic cell edge lengths: (1) 1.0, (2) 0.5 and (3) 0.25 mm 
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4 Comparison of Flat Plat Simulations 

Simulation Expense 

For the purpose of comparing simulations, the same simulation parameters listed in Table 2 are input 
to DAC and MONACO. Reference time-step and particle weight are determined by DAC’s preprocessor. 
These codes have different standard output processes. In DAC input file, the “number of time-steps 
between file updates” is set to 100 (setting it to 2000 resulted in 5.48 % increase in simulation real time). 
In MONACO input file, “Interval: Evaluate simulation process” and “Interval: Print output” are each set 
to 100, while “Interval: Write restart file” is set to 2000. 


Table 2. DSMC simulation parameters 


time-step 
particle weight 
transient period 
sampling period 
interval sample data 


3.102 X 10 -7 s 

4.645 x 10 9 real per simulated molecules 
8000 time-steps 
8000 time-steps 
2 
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Serial processing executables are generated with Intel R 10.1.018 compilers 34 using default 
optimization (-02) for serial processing. Each executable is run on an AMD Opteron™ Processor 244, 35 a 
64 bit processor having 1.8 GHz clock rate and 4 GB memory. Simulation time and memory are measured 
by Linux commands “time” and “top”, respectively. Table 3 summarizes simulation expense with serial 
processing. For this case, DAC utilizes 57 % less memory, while MONACO takes 9.8 % less time. 


Table 3. Serial processing simulation expense 


simulation code DAC MONACO 

real time (hr: min: s) 9:42:18 8:45:10 

memory 223 Mb 516 Mb 


Figure 4 demonstrates that the codes exhibit similar progression toward steady-state. Defining steady- 
state to begin when number of simulation molecules (simulators) remains consistently within two 
standard deviations from sampled population mean, then, DAC and MONACO simulator populations 
arrive at steady-state after 1500 time-steps, rounded to the nearest 100 th time-step. Number of gas-gas 
molecule collisions per time-step becomes steady after 1100 and 1000 time-steps, respectively. For both 
codes, number of gas-solid molecule collisions per time-step becomes steady after 700 time-steps. During 
sampling period, mean number of resident simulators in DAC simulation is 3,169,267 and 0.12 % 
greater than in MONACO simulation. Concurrently, mean number of gas-gas and gas-solid molecule 
collisions per time-step in DAC simulation are 50,042 and 894, respectively. These are 1.5 and 1.2 % 
greater than MONACO’S results, respectively. Simulations excluding molecule rotational energy 
exchange indicate that the minor differences in population and collision histories is mainly due to distinct 
molecule rotational energy exchange models. 
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Figure 4. Comparison of global simulator population histories 
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To further examine differences between code performance, they were recompiled with PGI R 10.4 
compilers 36 using default optimization (-01) and -Mprof=func for profiling of serial processing. A 
different compiler was selected for profiling because of profiling software availability. For profiling each 
code, the baseline simulation case is rerun under steady state for 20 time-steps with sampling. Figure 5 
compares temporal expense for the top eight expensive functions of each code. Percentages are relative to 
each code’s total run time. Simulation expense is dominated by molecule motion function “move” in both 
codes. Second most expensive function in DAC and MONACO is “read” and “calccells”, respectively. 
These functions indicate DAC is burdened by standard input; while, MONACO is burdened by grid 
managing processes. Third and fourth most expensive functions handle collision and random number 
generation processes, respectively. Other functions each take up less than 5% total run time. Overall, 
DAC’s function load balance is heavier on standard input and output processes and MONACO’S is 
heavier on grid managing processes. 
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Figure 5. Breakdown of simulation expense among most expensive functions of each code 



In addition, parallel processing performance is examined. For parallel processing, DAC and 
MONACO executables are generated with Intel® 10.1.018 compilers using default optimization (-02) 
with Open MPI 1.2. 7. 37 The selected simulation case is then repeated with 2, 4, 8, 16 and 32 processors, 
with each parallelized executable. Then, parallel processing speedup efficiency tj and memory spread 

S pm are computed from the following equations: 


r N t 

T proc 1, parallel 

T^roc 

S pm =— ( 2 ) 

m l 

where t\ and t para u e i are serial and parallel processing run time, respectively; N proc is the number of 
processors; m\ is the total processing memory of the serial run; and m pwc is the processing memory of 
processor proc of a parallel run. Serial processing time and memory are taken from Table 3. t pam u e i and 
m piV c are measured by “time” and “top”, respectively. Regarding parallel processing strategy, both codes 


assign one processor for standard input and output (I/O), however, DAC excludes this processor from 
DSMC work. 


Plots of i] ps and S pm constitute Figs. 6 (a) and (b), respectively. When N proc < 8, MONACO 

demonstrates greater speedup; however, when N proc > 8, DAC demonstrates greater speedup. This 
threshold is due to distinct prescribed master processor behavior and other less distinct differences in 
parallel processing code. Thus, speedup comparisons seem to favor the strategy of having master 
processor execute only standard I/O for N proc > 8. Regarding memory spread, MONACO demonstrates an 
erratic trend with increase in number of processors; whereas, DAC demonstrates a linear increase. Except 
case with 2 processors, memory spread resulting from DAC parallel processing is substantially less than 
that by MONACO. Greater memory spread by MONACO evidences inefficient parallel processing data 
handling procedures within that code. 




Figure 6. Comparison of parallel processing time and memory 


DSMC Metrics 

DSMC metrics gauge statistical validity. Generally, to have acceptable local simulator statistics, 
DSMC simulations abide by the following criteria: T / As > 1, r c /r 0 > 1, and N p /cell > 10. The 
present simulations obey these criteria throughout most of the computational domain. Small pockets 
ahead of leading edge and behind trailing edge are exceptions; however, because the study’s purpose is 
code comparisons, average grid converged surface properties are considered satisfactory. Figure 7 (a) 
compares field distributions of N v /cell. Qualitatively, both sets of contours are reasonable: minimum 
particle population occurs in the wake and maximum populations occur in compression regions. 
Deviations in population contours occur near horizontal wall surfaces away from leading and trailing 
edges, and in lower half of computational domain near the free stream. Simulations excluding molecule 
rotational energy exchange indicate that deviations near surface are due to distinct molecule rotational 
energy exchange models. However, deviations in lower half of domain, near free stream, also appear in 
those simulations. These deviations are due to asymmetry in MONACO results, which could be due to 
minor bias in selection of inflow and post collision velocities where y < 0. Determination of exact reason 
for this asymmetry is considered beyond the scope of this study. 
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(a) simulators per computational cell 
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Figure 7. Flow field distribution of DSMC metrics 


Figure 7 (b) compares field distributions of A/ As. Large differences in A/ As contours are due to 
differences in post-processing procedures. Although inconsequential for this case, it is relevant to note 
that DAC assumes As to be the average cell edge length, while MONACO assumes it to be the cubed root 
of the cell volume. These definitions are equivalent for cubes; however, they are distinct for general 
cuboids. Thus, in this case, the reason for distinct A/ As contours is due to distinct equations employed to 
compute A. While, DAC employs equation (1.12a) of Ref. 26, MONACO employs a form of equation 
(4.77) of the same reference. Because computational mesh is not adapted to A field, ascertaining the more 
accurate approach is considered outside the realm of this study. Nevertheless, it is notable that field 
distributions of A/As indicate that mesh adaptation according to MONACO’S result would require 
roughly twice the expense, throughout most of the domain, if cells are adapted to local mean-ffee-path. 

Table 4 lists the extrema of DSMC metrics. Large differences in values demonstrate sensitivity to 
differences in post-processing procedures. Extrema are associated with peak compression and expansion 
of flow in regions ahead of leading edge and behind trailing edge, respectively. Minimum mean-ffee-path 
A and collision time t c and maximum local molecule population N p are associated with compression 
region ahead of leading edge. Opposite extremes are associated with expansion region behind trailing 
edge. Extrema of DSMC metrics substantiate the merit of examining global field distributions of these 
metrics to evaluate overall simulation quality. Figs. 7 (a) and (b), indicate that DAC distribution of x c /x 0 
is reasonable. Adding procedures to extract this result in DAC is beyond the scope of this study. 
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Table 4. Extrema of DSMC metrics 


simulation code DAC 

X/As (min, max) (0.55,2.59) 

t c /t 0 (min, max) n/a 

N p /cell (min, max) (0,127) 


MONACO 

(0.59,432) 

(1.5,1513) 

(0.15,129) 


Flow Field Contours 

Figure 8 (a) compares contours of bulk flow velocity magnitude. Qualitatively, they indicate that the 
flow domain size is adequate. They demonstrate the presence of a boundary layer over the flat plate, 
located at 0 < x < 100 mm, and the presence of an oblique diffuse shock. Agreement in contours 
between the two codes is excellent, except near free stream, mostly, in lower half of domain. This 
deviation, again, reflects asymmetry in MONACO results. Similar results are observed with translational 
temperature contours. However, because of distinct molecule rotational energy exchange models, 
contours of rotational temperature demonstrate appreciable deviation throughout the entire domain, as 
illustrated in Fig. 8 (b). In summary, flow field contours demonstrate that MONACO is incorrect in the 
far field where y < 0. Both codes provide reasonable results in regions of significant gradients. 
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Figure 8. Comparison of (a) bulk velocity and (b) rotational temperature contours of N 2 flow field 
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Surface Properties 

Figures 9 (a) through (c) compare pressure, shear stress and heat flux along upper surface. Results of 
simulations performed without molecule rotational energy s rot exchange are also included. In addition, 
pressure and heat flux reported by Allegre and others 31 are included. Uncertainty of pressure and heat flux 
are +0.1Pa and +5 % q about expectation values. This uncertainty is based on measurement tools. 
Allegre and others did not report data relevant to uniformity of test section free stream flow. 

Surface pressure plots, Fig. 9 (a), indicate that differences in DAC and MONACO results are due to 
distinct s rot exchange models. Percentage changes are computed relative to DAC results. Mean 
percentage absolute-value-of change % | A | in surface pressure increases from 1.0 to 2.8 % as a 
consequence of including s rot exchange. For simulations that include s rot exchange, mean percentage 
absolute-value-of error %\Err\ of simulation results relative to laboratory expectation values is 12.4 and 
8.67 % for DAC and MONACO, respectively. Differences in accuracies are mostly due to DAC assuming 
constant rotational collision number, whereas, MONACO employs a variable rotational collision number 
model. 
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Figure 9. 


Comparison of upper surface properties 
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By contrast, for shear stress t x and heat flux q, the inclusion of s rot exchange decreases % | A | between 
the codes. Mean percentage change of t x decreases from 1.96 to 1.62 %, while that of q decreases from 
9.3 to 1.91 %. Thus, factors other than s rot exchange model cause differences between code simulation 
results. As indicated by field contour plots, these appear to be subtle differences in boundary and gaseous 
collision procedures. The unexpected percentage change in q between the codes for simulations excluding 
Srot exchange is suspected due to faulty values for s rot inserted by boundary condition or post processing 
procedures. Hence, q results, of simulations excluding s rot exchange, are not trusted. For simulations that 
include s rot exchange, % | Err | of simulation results relative to laboratory expectation values of q is 13.4 
and 12.8 % for DAC and MONACO, respectively. As expected, the variable rotational collision number 
procedure, on average, produces results closer to the laboratory expectation values. 

Flat plate drag, including pressure and shear components, are listed in Table 5. Given the 20:1 profile 
aspect ratio of the flat plate, it is reasonable to find drag due to shear 32 % greater than drag due to 
pressure. Percentage change in values between codes show that total surface drag is less sensitive than 
local surface aerothermodynamic properties by an order of magnitude. This substantiates employment of 
upper surface property distributions for comparing codes with selected simulation case. 


Table 5. Drag (N) of simulated flat plate 



D press 

D shear 

D 

DAC (DA) 

1.911xl0- 4 

2.472x10-4 

4.383x10-4 

MONACO (MO) 

1.916x10-4 

2.461x10-4 

4.377x10-4 

\MO — DA\/DA 

0.30 % 

0.46 % 

0.13 % 


5 Conclusion 

Two DSMC codes were compared: NASA’s standard production code DAC and an academic research 
code MONACO. These codes employ distinct grid processing procedures. They were compared by 
simulating a blunted flat plate wind tunnel test, using an identical flow field computational mesh. This 
mesh enabled satisfactory results according to a grid refinement study. 

Performance of simulations was compared. When performed on a single processor, DAC consumed 
57% less memory while MONACO consumed 9.8% less time. Function profiles of single processor 
simulations determined that DAC is characterized by more expensive standard input and output functions 
and MONACO is characterized by more expensive grid managing functions. Simulations were repeated 
with parallel processing with increasing number of processors. When number of processors was equal or 
greater than 8, DAC demonstrated greater speedup. This threshold is due to distinct prescribed master 
processor behavior and other less distinct differences in parallel processing code. Thus, speedup 
comparisons seem to favor strategy of having master processor execute only standard I/O for large 
number of processors; with threshold value expected to be case dependent. Regarding memory spread, 
MONACO demonstrated an erratic trend with increase in number of processors; whereas, DAC 
demonstrated a linear increase. Greater memory spread by MONACO evidences inefficient parallel 
processing data handling procedures within that code. In summary, parallel speedup and memory spread 
comparisons indicated that MONACO’S parallel processing procedures are generally less efficient than 
DAC’s. 
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Accuracy of simulations was also compared. DSMC metrics indicate that both codes yield acceptable 
quality solutions with same volume mesh and global weights. Flow field contours demonstrated that 
MONACO is incorrect in the far field where y < 0. Both codes provided reasonable results in regions of 
significant gradients. There, good agreement was observed in field contours of bulk velocity magnitude 
and translational temperature. Rotational temperature contours deviate throughout most of the domain 
because of distinct rotational energy exchange models. Namely, DAC assumes a constant rotational 
collision number, whereas, MONACO employs a variable rotational collision number model. For this 
reason, mean percentage error relative to laboratory expectation values of upper surface pressure and heat 
flux is slightly less for MONACO than for DAC. 

Overall, this study revealed that both codes, excluding grid adaptation, performed similarly. For 
parallel processing, DAC was generally more efficient. As expected, code accuracy was mainly 
dependent on the physical models employed. 
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